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Abstract 



This is a review of the problem of Mutually Unbiased Bases in finite di- 
mensional Hilbert spaces, real and complex. Also a geometric measure of 
"mubness" is introduced, and applied to some explicit calculations in six di- 
mensions (partly done by Bjorck and by Grassl). Although this does not yet 
solve any problem, some appealing structures emerge. 
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1. The problem 



A formula that has been of central importance in many discussions about 
the foundations of quantum mechanics is 



It expresses the complementarity of position and momentum. If we know 
everything about position, we know nothing of momentum. In Hilbert spaces 
of finite dimension N, the analogous equation would concern two orthonormal 
bases \e a ) and |/ ) such that 



The important thing is that the right hand side is independent of a and b. 
If such bases can be found, they are said to be mutually unbiased bases, or 
MUBs for short. The name emphasises that the information obtained in a 
projective measurement associated to one of the bases is completely unrelated 
to the information obtained from a projective measurement associated to the 
other. The question is: how many MUBs can one introduce in a given Hilbert 
space? If iV is a power of a prime number then one can find N + 1 MUBs 
[1] [2] , but for other values of N the number might be smaller. On the other 
hand it might not: I simply do not know. 

Why should you care about this problem? Apart from the fact that it 
is easy to state, several answers can be given. One answer is that MUBs 
are useful in quantum state tomography [3, 4]. Another has to do with 
various cryptographic protocols [5]. Thus, whether one wants to find or hide 
information, unbiasedness is a useful property. A third answer is that when 
one begins to look into it, the MUB problem leads one into many corners 
of mathematics that have been explored in communication theory, computer 
science, and so on, but which are relatively unknown to quantum physicists. 
If the essence of quantum mechanics is that it permits one to do things that 
cannot be done in a classical world, then many surprises may be lurking in 
those corners. 

As a matter of fact, over the past three years or so, papers about the 
MUB problem have appeared at a rate larger than once a month. I found 
about forty of those papers interesting, but I decided to quote only a small 



|(x|p)| 2 = constant . 
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fraction here. On the other hand I will borrow results freely from people that 
I mention in the acknowledgements; we recently wrote a paper containing 
much more detail than did my talk (details that this written version will 
occasionally hint at) [6]. 



2. Existing constructions 

Let us begin by taking a look at existing constructions. Recall that a discus- 
sion of position and momentum usually begins with the equation 

[x,p]=i, (3) 

or in terms of unitary operators 

In his 1928 book, Hermann Weyl [7] considered a finite dimensional analogue 
of this. We look for two unitary operators X and Z — the notation is supposed 
to suggest an analogy to the Pauli matrices — such that 

XZ = qZX , (5) 

where q is a phase factor. Weyl found that if q is a primitive root of unity, 
say if 



q = e 2 N , (6) 

then eq. (5) admits a representation that is unique up to unitary equivalence, 
and the eigenbases of X and Z are indeed mutually unbiased (although this 
piece of terminology came later!). If the eigenbasis of Z is chosen to be 
the standard basis, then the eigenbasis of X consists of the columns of the 
Fourier matrix F, whose matrix elements are 

F ab = q ab , < a, b < N - 1 . (7) 

It is called the Fourier matrix because it appears in the discrete Fourier trans- 
form. So far, all statements are independent of the dimension N. Closer 
examination of the group that is generated by X and Z reveals some dimen- 
sion dependent things; notably if N = p = a prime number, then one can 
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use the Weyl group to generate a set of N + 1 MUBs. But sometimes only 
three MUBs turn up in this way. The special status of prime numbers has to 
do with the fact that Weyl's representation theorem requires a phase factor 
that is a primitive root of unity, that is a root of unity such that q k ^ 1 for 
k < N. This is true for all the roots of unity only if N is prime. 

In general the following result has been established. Let N = p™ 1 -p^ 2 • . . . • 
pl k be the prime number decomposition of N, with p™ 1 < p 1 ^ < . . . < p^ k . 
Then the number of constructed MUBs obeys 

+ 1 < # MUBs < N +1 . (8) 

A complete set of N + 1 MUBs has been constructed for all prime power 
dimensions. Without going into any details (see Bandyopadhyay et al. [8] 
for a useful account), let me observe that the constructions differ somewhat 
depending on whether N is a power of an even or an odd prime. In par- 
ticular the behaviour of the MUBs under complex conjugation (relative to 
the standard basis) is strikingly different. I do not know what this means, if 
anything. 

When N is not a prime power, the known bounds are not very sharp. 
For some very special choices of N a somewhat higher lower bound is known 
[9], but in general the problem is wide open. For N = 6 several attempts to 
construct more than 3 MUBs have been made, using either group theoretical 
tricks or computer searches. There seems to be a growing consensus that 3 
is the best one can do. Perhaps this is related [10, 11] to the fact that no 
Graeco-Latin square of order 6 exists [12]. A set of N — 1 Latin squares that 
are mutually Graeco-Latin can be used to construct a finite afline plane. I 
leave this obscure remark as a hint that there are interesting connections 
between MUBs and combinatorics; for now let me just say that the available 
evidence concerning MUBs for A = 6 strikes me as rather weak. 

The MUB problem appears also in other branches of learning. Indeed 
the same, or nearly the same, problem has occurred in the theory of radar 
signals [13], to operator algebra theorists [14], in Lie algebra theory [15], and 
in coding theory [16]. Unfortunately all technical terms differ between these 
groups of authors, so it is quite difficult for a quantum theorist to extract 
the information gathered in any of the other fields. Suffice it to say that the 
Lie algebra theorists have the best name for the problem. They call it the 
"Winnie-the-Pooh problem" , for reasons that would take us too far afield. 
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(It has to do with a free Russian translation of a verse hummed by Pooh.) 
The actual results achieved in the various fields appear to be roughly the 
same, as far as I can tell. 



3. A packing problem 

How should we look at the MUB problem? A first way is to view it as a 
packing problem in Hilbert space, or more accurately in complex projective 
space. 

In Hilbert space a basis can be represented as the columns of a unitary 
matrix. Assuming one basis is represented by the unit matrix, all bases that 
are MUB with respect to it must be represented by unitary matrices of the 
form (for N = 3, say) 



U = 



1 


1 


1 


g*01O 


g*011 


g*</>i2 


g#20 


gi'fel 


g#22 



(9) 



This is (except for the normalizing factor) a complex Hadamard matrix. 
The first row has been choosen to contain only ones by convention. All the 
vectors in all bases that are MUB with respect to the standard basis can 
therefore be found on a torus parametrized by iV — 1 phases. This torus 
has a natural interpretation as the maximal flat torus in complex projective 
space (equipped with the Fubini-Study metric). Now I seem to be saying 
that finding the MUBs is equivalent to a packing problem on a flat torus, but 
unfortunately this is not quite true, because the tori in complex projective 
space are not totally geodesic. What this means is that intrinsic distance 
on the torus does not directly reflect the Fubini-Study distance [17]. We do 
have a packing problem in complex projective space, but packing problems 
are difficult, and moreover their solutions tend to depend on dimension in 
peculiar ways. 

A straightforward approach, while we remain in the N dimensional Hilbert 
space, is to begin by asking for a classification of all complex Hadamard ma- 
trices. But here the existing results are very incomplete [18, 19, 20]. For any 
N, the Fourier matrix exists, hence one basis that is MUB relative to the 
standard basis always exists. Are there more? We can multiply the Fourier 
matrix from the left with a diagonal unitary matrix, that is to say we can 
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multiply the rows with suitable phase factors, and hope that with appro- 
priate choices of the phase factors we can find further bases that are MUB 
relative to both the standard basis and to the Fourier basis. We can also ask 
if there are further Hadamard matrices, not related to the Fourier matrix in 
this way. One would naively expect that the MUB problem should become 
easier the more such Hadamards one finds. This expectation is not borne 
out however. 

To classify (complex) Hadamard matrices, it is helpful to begin by declar- 
ing two Hadamard matrices to be equivalent if one can be reached from the 
other by permutations of rows and columns, and by multiplying rows and 
columns with arbitrary phase factors. Since we think of the columns of an 
Hadamard matrix as representing a basis in Hilbert space, these operations 
applied to columns mean nothing at all to us, while the same operations 
applied to rows are achieved by unitary transformations that leave the stan- 
dard basis invariant. I will refer to the standard basis as the "zeroth MUB". 
Then I do not loose any generality if I assume that the first vector in the 
"first MUB" has entries l/y/N only. Thus (MUB)! is represented by an 
Hadamard matrix in dephased form (the first row and the first column has 
all entries real and positive), while any further MUBs are represented by 
enphased Hadamard matrices. 

When N = 2 or 3 the Fourier matrix is unique up to equivalences. Hence, 
once (MUB) and the first vector of (MUB)i are chosen, everything else is 
forced. But when N = 4 there exists a one parameter family of inequivalent 
Hadamard matrices, found (appropriately) by Hadamard himself. It is 



" 1 1 


1 


1 " 


1 e** 


-1 


_ e i<t> 


1 -1 


1 


-1 


. 1 -e^ 


-1 


e i<t> 



(10) 



Hence there is some freedom in choosing (MUB)!. And the choice matters — 
we have to set (j) equal to zero, or to n, if we want to have an additional three 
MUBs. 

When N = 5 the Fourier matrix is again unique [18], but when N = 6 
there are many choices. I will discuss them later. When N is a prime number 
one cannot introduce any free parameters into the Fourier matrix, but at least 
when N = 7 another, unrelated, one parameter family is known. In brief, 
the situation is confusing, but based on the information available one would 
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be inclined to guess that finding MUBs should be particularly easy when 
N — 6. Which is definitely not the case! 



4. The shape of the body of density matrices 

A second way to look at MUBs is to look at density matrices. Each vector in 
an N dimensional Hilbert space can be thought of as a rank one projector in 
the set of Hermitian matrices of unit trace, which has real dimension N 2 — 1. 
It is convenient to think of this space as a vector space, with the origin placed 
at 



Then a unit vector |e) in Hilbert space corresponds to a real unit vector e in 
an N 2 — 1 dimensional vector space, whose elements are traceless matrices. 
The explicit correspondence is 



The catch is that only a small subset of all unit vectors in the large vector 
space arises in this way. Anyway, the body of density matrices is now ob- 
tained by taking the convex cover of all those points on the unit sphere where 
some vector e ends. Equivalently, the body of density matrices consists of 
all Hermitian matrices of unit trace and positive spectrum. Either way, it 
is a body with an intricate shape that is difficult to visualize. It touches its 
N 2 — 2 dimensional outsphere in a small 2(iV — 1) dimensional subset, arising 
from vectors in Hilbert space through the above correspondence. These are 
the pure states. The whole space is naturally an Euclidean space, with a 
squared distance given by 




(12) 



Our vector space is equipped with the scalar product 



ef 



-Tref . 

2 



(13) 




(14) 
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This is the Hilbert-Schmidt distance. The geometry induced on the set of 
pure states is of course precisely the Fubini-Study geometry on complex pro- 
jective (N — l)-space. Alternatively, the Hilbert-Schmidt distance provides 
the chordal distance between two points on the outsphere of the body of 
density matrices [17]. 

An orthonormal basis in Hilbert space is now represented as a regular sim- 
plex with N corners, spanning an (N — 1) dimensional flat subspace through 
the center of the outsphere. Two bases will be MUB if and only if their 
respective (N — l)-planes are totally orthogonal. Since the whole space has 
dimension 

N 2 -l = (N+1)(N-1) , (15) 

it is clear that at most N + 1 MUBs can be found. 

After a moment's reflection one sees that a complete set of MUBs defines 
a rather interesting convex polytope, which is the convex cover of N + 1 
regular simplices placed in totally orthogonal (N — l)-planes. It is called 
the Complementarity Polytope [21]. Evidently such a polytope will exist in 
every iV 2 — 1 dimensional flat space. This does not yet solve our problem 
though, because given such a polytope its corners will correspond to vectors 
in Hilbert space if and only if we can rotate it so that all its corners fit into 
a special 2(N — 1) dimensional subset of its N 2 — 2 dimensional outsphere. 
And when N > 2 this is hard! 

Nevertheless this is an appealing picture of MUBs. In particular, this is 
the way to see why MUBs solve the problem of optimal state determination 
in non- adaptive quantum state tomography [1, 2]. 

5. Real Hilbert spaces 

It is instructive to pause to think about MUBs in real Hilbert spaces, because 
in this case it is easy to derive some negative results. Indeed real Hadamard 
matrices can exist only if N is two or divisible by four. So we see at once 
that in a three dimensional real Hilbert space, it is impossible to find even 
a pair of MUBs. We can see why geometrically — in real Hilbert space, not 
among the density matrices. If a vector is represented by a pair of antipodal 
points on the unit sphere in R 3 , then a basis is represented by the corners of 
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an octahedron. It is then geometrically evident that there are four vectors 
that are unbiased with respect to a given basis, and they form the corners of 
a cube that is dual to the octahedron. But they do not form a basis! 

In R 4 the situation is different. It is still true that a basis is represented 
by the corners of a cross polytope (the generalisation to arbitrary dimension 
of the octahedron), and there will be eight vectors that are MUB with respect 
to a given basis, again forming a cube that is dual with respect to the cross 
polytope. But specifically in four dimensions, a cube can be regarded as the 
convex cover of two symmetrically placed cross polytopes [22]. In this way 
we end up with three bases represented by three symmetrically placed cross 
polytopes, and their convex cover is a famous Platonic body known as the 
24-cell (having no analogue in three dimensions). This gives us three MUBs, 
and since the dimension of the set of real four-by-four density matrices is 
nine, this is a complete set in four dimensions. It is a set that has acquired 
some fame in quantum foundations, because a pair of dual 24-cells correspond 
to 24 vectors that are uncolourable in the Kochen-Specker sense [23]. 

If you want to know what happens in higher real dimensions, consult 
Boykin et al. [24]. 

6. Mubness 

Let us return to the complex case. I will describe an attempt to investi- 
gate what happens in six dimensions, but in order to say something more 
interesting than the obvious "I failed" , I need a measure of how much I fail. 
Many such measures, of varying degrees of sophistication, can be imagined. 
The one we use is based on the picture of MUBs that emerged from the den- 
sity matrix point of view, namely that they correspond to totally orthogonal 
(N — l)-planes in an N 2 — 1 dimensional space. Just as vectors in an N 
dimensional space can be regarded as rank one projectors in a higher dimen- 
sional space, so one can regard n-planes in an m-dimensional space as rank n 
projectors in a vector space of sufficiently high dimension. In mathematical 
terms, this provides an embedding of the Grassmannian of n-planes into a 
flat vector space. The rank n projectors will sit on a sphere in this flat space, 
and its natural Euclidean distance provides us with a chordal distance be- 
tween the projectors. This notion of distance has been used to study packing 
problems for n-planes [25], and it is the one we use to measure the distance 
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between bases in Hilbert space. The chordal distance attains its maximum 
if the n-planes are totally orthogonal, that is to say, if the bases are MUB. 
In this way it does provide a measure of "mubness" . 

The details, adapted to our case, are as follows. Starting from a basis in 
Hilbert space, form the N vectors e a . Then form the (iV 2 — 1) x N matrix 



B = 



IN 



N 



■[ei e 2 ... e N ] . 



(16) 



It has rank N — 1. Next we introduce an (N 2 — 1) x (N 2 — 1) matrix of fixed 
trace, which is a projector onto the (JV — 1) dimensional plane spanned by 
the e„: 



p = B B l = 



N- 1 



N 



[d ... e N ] 



L C N J 



(17) 



It is easy to check (through acting on e a say) that this really is a projector. 
The chordal Grassmannian distance between two N — 1 planes spanned by 
two different bases then becomes 



D 2 (P 1 ,P 2 ) = iTr(P 1 



P 2 ) 2 = N - 1 - TrP x P 2 . 



(18) 



It should be obvious that there is an analogy to how the density matrices 
were defined in the first place, and to the Hilbert-Schmidt distance between 
them. The difference is that now the projectors represent entire bases in 
Hilbert space, not single vectors. 

Working through the details, one finds that 



P 2 (P!,P 2 ) — N — 1 - EE (l(^IA)l 2 - 



(19) 



Thus 



< D z c < N - 1 



(20) 



and the distance attains its maximum value if and only if the bases are 
MUB. As a measure of mubness, the chordal Grassmannian distance has 
the advantages that it is geometrically natural and simple to compute. It is 
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also natural from the tomographic point of view, although I certainly cannot 
claim any precise operational meaning for it. 



7. N = 6 

My third and final way to look at MUBs is to simply perform calculations to 
see what happens, without thinking very much. We tried the first open case: 
N — 6. Maybe this was a mistake, because 6 clearly sits astride the even and 
the odd prime numbers. Perhaps we should concentrate on 15 = 3 • 5? (The 
question how many mutually Graeco-Latin squares exist is actually open in 
this case.) On the other hand, a six dimensional Hilbert space is already a 
very large space to search in — and fifteen is larger. So we stick to six. 

Let (MUB)o be the standard basis, and (MUB)i be given by the columns 
of some dephased Hadamard matrix (that is one whose first row and first 
column are real). We make a choice for (MUB)i, find all enphased Hadamard 
matrices that represent bases that are MUB with respect to (MUB)i, and 
then check how far apart the latter are, in the sense of the chordal distance. 

When N = 6 there are several choices for (MUB)i. The following de- 
phased Hadamard matrices are known [20]: 

i) The Fourier matrix, augmented with two free parameters: 



" 1 


1 


1 


1 


1 


1 


1 




gV* 2 


g 3 


1 4 e i0i 


q 5 e i<h 


1 


q 2 


q 4 


1 


q 2 


q 4 


1 




g«02 


g 3 




g3gi</>2 


1 


q 4 


q 2 


1 


q 4 


q 2 


_ 1 




q 4 e i<t> 2 


g 3 


q 2 e l<j)1 


qe t<t>2 



There are various equivalences of the form F(0!,0 2 ) ~ -^(03?04); this has 
been sorted out, but this is not the place to give all the details. 

ii) The transpose of the above, again with two free parameters. 

iii) A one-parameter family of matrices D(<f>) whose entries, when the free 
phase is set to zero, are fourth roots of unity. It is known as the Di^a family. 

iv) An isolated matrix S whose entries are third roots of unity. 
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v) A matrix found by Bjorck, which is most conveniently given as the circu- 
lant matrix 



C 



1 


id 


-d 


—i 


—d id 


id 


1 


id 


-d 


—i — d 


-d 


id 


1 


id 


—d —i 


—i 


-d 


id 


1 


id —d 


-d 


—i 


-d 


id 


1 id 


id 


-d 


—i 


-d 


id 1 



(22) 



A matrix is said to be circulant if its columns are cyclic permutations of its 
first column. The complex number d has modulus unity and is 



d = 




^ => d 1 



[1 - V3)d+ 1 = . 



(23) 



Several people have invested some effort in making this list as long as it can 
be, and up to two weeks before my talk, I thought that it might well be a 
complete list. Still, we do have considerable latitude in how we choose the 
first MUB. 

Let us begin with 



(MUB)i = F(0, 0) . (24) 

It happens that all vectors that are unbiased with respect to this choice of 
the zeroth and first MUB have been computed, first by Bjorck and Froberg 
[26], and independently by Grassl [27]. Bjorck did not express the problem 
in these terms however. He was interested in biunimodular sequences, that is 
to say sequences of unimodular complex numbers x a whose discrete Fourier 
transform 

^ = ^=Y^x b q ba (25) 

also consists of unimodular complex numbers. On reflection, one sees that 
the two problems are equivalent. Bjorck and coauthors eventually solved this 
problem for all N < 8. When = 6 there are altogether 48 such sequences; 
12 Gaussian ones — they were known to Gauss — and an additional 36. The 
Gaussian ones have entries that are 12th roots of unity, while the additional 
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ones involve Bjorck's magical number d. See also Haagerup [18], who seems 
to have been the first to get this quite right. 

More is true. For a biunimodular sequence the autocorrelation function 

is 

j N-1 y N ~ 1 

7b = T7 Yl XaXa+b = H l^al 2 ^"^ = $6,0 • (26) 
iV a=0 iV a=0 

Therefore x a and x a +b, with b fixed and non-zero, are orthogonal vectors. 
Then it follows that the 12 + 36 vectors found by Bjorck can be assembled 
into 2 + 6 unitary circulant matrices that are MUB with respect to the 
standard basis and the Fourier basis. When Grassl redid this calculation 
(using the program MAGMA) he observed that each of the 48 vectors can 
be used to form a basis in exactly two ways. Thus we end up with exactly 2 
+ 2 + 6 + 6 = 16 possible choices for (MUB) 2 . In itself, this is more than 
we need for a complete set of 7 MUBs. Provided that the Fourier basis is 
included, the question whether one can find a fourth MUB boils down to 
the question whether the chordal distance squared between any pair among 
the 16 is equal to 5, the maximal distance squared attained by MUBs in six 
dimensions. 

The answer is no. The 4 Gaussian MUBs, composed of 12th roots, form 
a perfect square with side lengths squared D\ = 2. The two groups of 
"non-Gaussian" MUBs are isometric copies of each other. One group con- 
sists of circulant matrices, while the other group consists of Fourier matrices 
enphased using the magical number d. The distance squared between any 
Gaussian and any non-Gaussian MUB is D 2 C « 4.62 — rather close to 5, if one 
takes an optimistic view of things. The distance between the two six-plets of 
non-Gaussian MUBs is always D 2 « 3.71. Inside each group, distances reach 
all the way up to D 2 ps 4.64, which is even closer to 5. 

Still, although the pattern is nice, the conclusion is negative: the Fourier 
matrix cannot be included in a set of more than three MUBs. We do not 
have complete results for other choices of the first MUB. We made a program 
that lists all MUB triplets where all the entries of the matrices are 24th roots 
of unity. Quite a few triplets, with quite interesting structures, did turn up 
in this way, but there were no MUB quartets. Note that in prime power 
dimensions the standard solution for complete sets of MUBs contain iVth 
or 2iVth roots of unity only (depending on whether N is odd or even [2]); 
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in other words the analogous calculation in arbitrary dimension would have 
found the known complete sets. 

Although the evidence is incomplete, we do seem to be driven to admit 
that there can exist at most 3 MUBs when N = 6. A nagging doubt remains, 
because there is always the possibility that the above list of Hadamard matri- 
ces is incomplete. In particular, could the parameter spaces be incomplete? 
We do know that the number of free parameters in Fourier's, Diva's, and 
Bjorck's matrices is at most 4, while the matrix S cannot have any free 
parameters at all [28]. But we do not know if there are that many free 
parameters. 

While I was thinking about what to say in my talk, Beauchamp and 
Nicoara [29] intervened. They found a new one-parameter family of Hadamard 
matrices, namely 
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-1 


_ i 
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-y 
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—x 
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-l 
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1 
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1 
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-t 
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_1 




y 
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X 


-t 


t 


—x 


-1 



B 



where (x, y, z, t) are complex numbers of modulus one, related by 



(27) 



t = xyz 



(28) 



x 



l + 2y-y 2 
y(-l + 2y + y 2 ) 

l + 2y + y 2 ± V2y/1 + 2y + 2y i + 



(29) 



(30) 



1 + 2y - y 2 

with y remaining as a free phase factor. By construction this family contains 
all iV = 6 Hermitian Hadamard matrices. The phase of y cannot be chosen 
quite arbitrarily; an interval around y — 1 is excluded. On closer inspection 
one finds that this family starts from a matrix that is equivalent to Bjorck's, 
passes through the Di^a family, comes back to Bjorck, repeats twice, and ends 
at Bjorck. The two branches of the square root lead to equivalent families. 
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What does this mean? I do not know. If it means that there are four 
dimensional families of Hadamard matrices, including the Fourier matrix and 
Bjorck's matrix, then it also means that we have looked in a very small part 
of parameter space only. In fact we do have some reasons to believe that this 
is really so. Therefore it seems to me that the conclusion that we must draw 
about the MUB problem in six dimensions is: We have almost no evidence 
either way. 

8. The real problem 

The distance that we introduced can, in principle, be used to convert the 
MUB problem into that of maximising a function, such as 

F = Y,D 2 c (P l ,P J ). (31) 

A similar procedure has been used [30] to find approximations to SIC-POVMs — 
this particular acronym stands for a kind of relative of the MUB problem — in 
dimensions up to iV = 45. In our case the upper bound is attained if the 
N + 1 projectors represent totally orthogonal (N — l)-planes. Whether we 
can reach the upper bound using (N — l)-planes spanned by bases in the 
underlying Hilbert space is of course precisely the question. 

I should add that I have not really done justice to the point of view that I 
tried to stress in the beginning, that the MUB problem leads one into many 
corners of useful mathematics that have not been very much explored by 
quantum physicists. But if you search your favourite eprint archive for some 
of the many papers, whose existence I hinted at, you will see what I mean. 
Meanwhile, the somewhat botanical spirit of my talk is perhaps appropriate 
in the town where Linnaeus was educated. 

Anyway the real "MUB problem" is not how many MUBs we can find. 
The real MUB problem is to find out what we can do with those that exist. 
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